The Tomonaga-Luttinger Model and the Chern-Simons Theory for the Edges of 
Multi-layer Fractional Quantum Hall Systems 
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Wen's chiral Tomonaga-Luttinger model for the edge of an m-layer quantum Hall system of total 
filling factor v = m/(pm ± 1) with even p, is derived as a random-phase approximation of the 
Chern-Simons theory for these states. The theory allows for a description of edges both in and out 
of equilibrium, including their collective excitation spectrum and the tunneling exponent into the 
edge. While the tunneling exponent is insensitive to the details of a v — m/(pm + 1) edge, it tends 
to decrease when a v = m/(pm — 1) edge is taken out of equilibrium. The applicability of the theory 
to fractional quantum Hall states in a single layer is discussed. 

PACS numbers: 71.10.Pm 
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I. INTRODUCTION 

The physics at the edges of quantum Hall systems has 
been the subject of both theoretical and experimental 
research since the early days following the discovery of 
the integer quantum Hall effect*, and up to the present. 
When the system exhibits the integer or the fractional 
quantum Hall effect (FQHE) the edge attains a unique 
role as the only region in the sample where gapless excita- 
tions exist. The low-energy physics associated with these 
edge excitations is especially rich under the conditions of 
the FQHE. In a pioneering series of papers^ Wen has de- 
veloped a theory for the edges of FQHE samples in terms 
of several interacting chiral Tomonaga-Luttinger liquids 
(CTLL), and predicted a power-law current- voltage char- 
acteristic for tunneling into such edges. 

Subsequent tunneling experiments^* 6 * 7 - have found a 
non-ohmic conductivity I ~ V 1 , with an exponent which 
approximately (and non universally 7 ) scales as 7 ~ \/v. 
However, with the exception of filling factor v = 1/3, this 
result is at odds with the predictions of Wen's theory. 
Moreover, the experiments suggest that a single charged- 
mode CTLL model appears to apply to filling factors in 
which the system is compressible and lacks the energy- 
gap needed to distinguish the edge as the region where 
the low-energy physics resides. While various theories 
have been put forward there is no consensus on the solu- 
tion to the problem 8 . 

In order to better understand the conditions under 
which the CTLL model for the FQHE edge is appli- 
cable it is desirable to test its predictions and study 
the way it arises from the underlying theory using as 
many theoretical approaches as possible. Previous re- 
search aiming at achieving this objective has included 
numerical calculations2*i2*ii as well as analytical studies 
with Laughlin's wave- function^ and the Chern-Simons 
(CS) theorjiisliiii for the FQHE as their departure points. 
These studies have been concerned predominantly with 
the Laughlin fractions v = l/(p + 1), with even p. Here 
we wish to extend our previous treatment of Ref. ^3 
to multi-component FQHE systems. We concentrate on 
m-layer quantum Hall systems with total filling factor 



v = m/(pm ±1), with even p, whose state is described 
by an m-component generalization of the Laughlin wave- 
function, as first introduced by Halperin 1 ^. We also con- 
sider the issues which arise when trying to apply a similar 
treatment to the Jain fractions, i.e. to the single-layer 
FQHE states with the same filling factors. 

Our plan is to obtain the one-dimensional CTLL the- 
ory for the edge as the low-energy limit of the underly- 
ing CS theory for the two-dimensional system. The CS 
description, presented in Section [H] follows closely the 
points of view taken by Jain 16 and Wenii. We first trans- 
form the electrons to composite fermions which fill, for 
the special filling factors under consideration, the lowest 
effective Landau level in each layer of the multi-layer sys- 
tem or m such levels in the single-layer problem. We then 
transform each of the layers (or in the case of the Jain 
states, filled Landau levels) into a bosonic condensate. 
The study of the coupled dynamics of these condensates 
constitutes the major part of the work. 

We begin, in Section lHll by considering the solutions to 
the CS mean- field equations. We find that the solutions 
form an m-dimensional manifold parameterized by the 
conserved momenta along the edge of the m condensates, 
or cquivalently, by the charges which they carry. Only a 
small fraction of the manifold describes equilibrium con- 
figurations which minimize the energy for a given total 
charge at the edge. These configurations are typically, 
but not always, symmetrical with all the condensates 
sharing the same state. The remaining solutions cor- 
respond to edges which are not in equilibrium and are 
generally asymmetrical. The conditions for the realiza- 
tion of a particular solution are discussed. 

Once the mean-field configurations are at hand we pro- 
ceed to consider small oscillations about them and obtain 
the edge excitations as the random-phase approximation 
(RPA) modes of the theory. We identify m RPA branches 
of edge excitations whose k — > limit is composed of the 
edge states of the m types of vortex-like excitations that 
exist in the bulk. Out of these m modes one is the edge 
magnetoplasmon which arises from in-phase oscillations 
of the condensates while the remaining posses acoustic 
dispersion, propagate with or against the magnetoplas- 
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mon [for v = m/(pm + 1) and v — m/(pm — I) respec- 
tively] and describe out of phase motion of the electron 
liquids. While most characteristics of the edge modes 
do not depend on the detailed nature of the underlying 
mean-field solution there are, however, some differences 
between the excitations of symmetrical and asymmetri- 
cal configurations; for example the net charge which is 
involved in the fluctuations. Consequently, we consider 
the two cases separately in Sections |^ and 

Finally, by using the edge modes to expand the devi- 
ations of the various fields from their mean-field average 
configurations we recover the CTLL model for the edge. 
The expansion coefficients become upon quantization the 
density operators in terms of which the model is formu- 
lated. Their commutation algebra as well as their Hamil- 
tonian are deduced from the quadratic expansion of the 
CS theory in the deviations. By inverting the CS trans- 
formation we also obtain the bosonization formula which 
expresses the original electronic operator at the edge in 
terms of the density operators. These results can be used 
to calculate the tunneling exponent into the edge. In ac- 
cord with the results of Weni£, we find that the exponent 
is insensitive to the details of the edge configuration or 
the value of m when v = m/(pm +1). On the other 
hand, it tends to decrease when a v — m/(pm — 1) edge 
is taken out of equilibrium. 



II. THE CHERN-SIMONS THEORY 

In the following we consider a semi-infinite sample sub- 
jected to a constant magnetic field Bz, with an edge de- 
fined by an external potential Aq(x) and periodic bound- 
ary conditions over length L in the y direction. The elec- 
tronic spin is assumed to be polarized by the magnetic 
field and will be left out of the analysis. 

We are interested in studying a particular correlated 
state of the m-layer system in which each of the layers is 
in a v — l/(pm ± 1) state with strong inter-layer corre- 
lations. An appropriate wave-function to describe such 
a system was introduced by Halperin 15 . For a two-layer 
system it is of the form 



* oc J] ( Zla - z l0 y ±l Yl (z 2a - z 20 y ±x 

a<(3 a<{3 



(1) 



where Zi a is the complex coordinate (in units of the mag- 
netic length I = y/hc/eB) of electron a in layer I. 

An alternative field theoretical description of the same 
states may be achieved via the CS transformation from 
the original electronic operators on the m layers, to 
new operators, ipi, describing composite fermions made 
out of p flux quanta (j>Q — he/ e attached to each electron^ 



^( r ) = e 'AW^( r ) , A(r) = / dV/i(r - r')p(r') 



(2) 

where p(r) = */( r )*/( r ) = J2i ^/( r )^/( r ) is tne 
total density. The new fields are fermionic provided 
/i(r — r') = /i(r' — r) +pn and p is even. The CS trans- 
formation introduces additional terms in the equation of 
motion of the which couple to them as vector and 

scalar potentials 



a(r) = ~ ZdVF^r-rOpCr'), (3) 
a (r) =-~J dVF x (r - r>(r>(r') , (4) 

where Fi(r — r') = V, /i(r — r') and v(r) is the total ve- 
locity operator. These " statistical gauge potentials" tend 
to reduce the external fields. As a result, for electronic 
filling factors v = m/(pm± 1) and within the mean-field 
approximation, the composite fermions filling factor in 
the residual magnetic field is vcf — ± m with a single 
occupied Landau level in each layer [the minus sign cor- 
responds to the fractions v — m/ (pm — 1) for which the 
direction of the residual field is opposite to that of the 
external field]. 

The resulting Lagrangian density reads^S^ 



1 

2M 

e 



-iHV + - (A + a) 
c 



0/ 



-\ f dV P (r)[/(r - r')p(r') - , 
2 J 2p(f>o 



(5) 



where M is the band mass. Here we have assumed that 
there is no tunneling between different layers. We will 
comment on the effects of such tunneling processes when 
we discuss the electronic propogator at the edge. 



For a single-layer v = m/(pm ± I) FQHE sys- 
tem the above transformation results, on the mean- 
field level, in ±m filled composite fermions Landau lev- 
els within the layer. One can use the Landau wave- 
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functions in the residual field (and in the presence of 
the edge) to expand the composite fermion operator 
V^(r) = J2 Iik c Iik L-V 2 e ik y Ul (x - kl*) = £ f tfj(r). 
The operators ipj defined here are not fcrmionic. Al- 
though operators which correspond to different Landau 
levels anti-commute, operators on the same level obey 
W>/(r),^(r')} = SuA^x^y-y') ± 5 TJ 5(r - r'). 
For example, away from the edge Ai(x,x\y — y') — 



l/(27rZ' 2 )exp[-(l/4)(r- 



-(i/2)(x+x')(y—y')], where 



r is measured in units of the magnetic length in the resid- 
ual field I' = y/m/vl. Since the extent of A/ in the bulk 
is roughly \fll' and as long as the length scales under 
consideration are larger than V and I is not too large, 
one may try to approximate the anti-commutator (after 
an appropriate rescaling of the fields) by a delta-function. 
It should be noted, however, that the integral over r of 
Aj vanishes for even I and it is not clear to what extent 
such an approximation is valid (see also in this context 
Ref. l20Th Ignoring this issue one is able to recover the 
Lagrangian density © provided that off-diagonal terms 
in the index /, corresponding to transitions between dif- 
ferent Landau levels, are neglected. Since these terms are 
not smaller than the diagonal ones such a step is not a 
priori justified. 
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In order to simplify the analysis we will assume that 
the confining potential is sharp enough so that the filled 
Landau levels extend from the edge all the way into the 
bulk (see Fig. Cleaved edge overgrown samples ap- 
pear to be good realizations of this limit 2 i. As the con- 
fining potential is made smoother higher Landau levels 
may cross the Fermi energy in the vicinity of the edg o 22 i 23 
or the edge may reconstruct by creating a strip of elec- 
trons separated from the bulk of the sampl o 24 ! 25 . Under 
such conditions one expects to find extra edge channels 
and edge modes. These cases will not be addressed in 
the present work and we will consider m electronic liq- 
uids (one per layer) that occupy the entire plane of the 
system. 

Next, we transform each of the fermionic fields ipi 
a bosonic field <f>i through an attachment of a single flux 
line in the ±z direction, namely, 

«Mr) = e iA '«^(r) , A,(r) = / d 2 r> / 2 (r - r>/(r') , 



o.o 



2 

-kf 



(6) 

with /a(r — r') = / 2 (r' — r) ± tt. The new gauge poten- 
tials introduced by this transformation 

aj(r) = / dVF 2 (r - r') Pl (r') , (7) 
a Jo (r) = T-J dVF 2 (r - r') Pl (r')vj(r') , (8) 

where F 2 (r — r') = V r / 2 (r — r'), cancel on the average 
the residual magnetic field in the bulk of the sample. 
Their dynamics is controlled by an additional CS term 
=F(e/2(/> ) Ej=i e ^ l/c ' a i„dv a i a m the Lagrangian density 
JSJ). Note that contrary to the flux lines attached by 
the first transformation J5J the new ones are felt only by 
particles in the same condensate as the particle carrying 
the flux. 

By shifting to eij + a M , integrating over and 
introducing the phase and velocity fields 



FIG. 1: A schematic plot of the composite fermions Landau 
levels for a sharp v — 2/3 edge. The dotted line represents the 
Fermi energy and k denotes the conserved y momentum of the 
states. When the edge is made smoother higher Landau levels 
(indicated by the dashed line) may cross the Fermi energy at 
the vicinity of the edge. 



6z = Vp7e^, v T = Avfr+^aj + A), (9) 



we obtain the final hydrodynamic form of the Lagrangian 
and Hamiltonian densities 26 



i=i 



PI (Ac + a Io - d t 6j) - \ Pl v] l -YPlL 



i m 1 r 

^ £ XuPiPj-^ J dVp(r)[/(r-r>(r') 



I,J=1 



m r - 

£ (K- 1 ) IJ -e^dtv] (2d 3 9j - 4) + a Io (e^ - l) 



j,j=i 



(10) 
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1=1 L 



PiA Q 



i - i r 

^ £ WW + ^ J d\' P {v)U{r 



i,j=i 

r 



r')p(r') 



(11) 



Henceforth length is measured in units of the magnetic 
length I, time is normalized by the inverse cyclotron fre- 
quency uj c = eB/Mc, energy by hu> c and the density by 
the Landau level degeneracy po = l/2nl 2 . The above 
densities have been derived under the assumption that 
no localized vortices are present in the sample. An- 
ticipating the use of the mean-field approximation we 
have included in (jlOll ip> additional contact interactions 
in order to partially account for the kinetic energy com- 
ing from the short range structure of the wave-function 
which is absent in this approximation^ 7 *^. We com- 
ment on the values of their couplings Ajj below. The 
matrix (if -1 )^ = ^fp/(pm ± 1) ± §u is the inverse of 
K = pC ± 1 with 1 and C being the identity and pseudo- 
identity matrices (Cu = 1). The matrix K is an example 
of the general classification scheme, introduced by Wen 
and Zeo 17 i 30 , of Abelian FQHE states. Its elements Ku 
are the number of flux lines affecting the phase of a par- 
ticle in the 7-th condensate when it encircles a particle 
in the J-th condensate. 



III. MEAN-FIELD SOLUTIONS 

Variation of the action l|l(J|) with respect to the vari- 
ables ai ,vi,9i and pj results in the mean-field equations 



£(# l ) IJ e ij d i v 3 J = ^- Pl , 



(12a) 



j=i 



£ {K- l ) u {dtdiOj - d t v]j - d iajo ) = e ijPl vi , (12b) 
J=i 

d t pi = -V(piui) , (12c) 

d t 6 I = - 1 -vj+ 1 -^l+A + a Io 
2 2 ^Jpi 

-^W/-^ / t/(r-r>(r')dV.(12d) 

In the absence of confining potentials these equations 
support a uniform solution pi = v/m with Vj = <Zj = 
and 6j = —pt. The correct noninteracting (U = 0) en- 
ergy density Tt = v?ku c /2 for this solution is recovered by 
taking Xu = 2irm/ v5u . One also finds pi = hw c indi- 
cating that mean-field configurations with slightly differ- 
ent filling factor and which are continuously connected to 
the constant solution necessarily involve higher Landau 
levels^. This choice is not unique. Taking Ajj = 2ttKjj 
also recovers the correct noninteracting energetics. As 
we shall see, the form of Xjj affects some of the details 
of the mean field solutions. It does not, however, modify 
any of the general conclusions of our study. 



Up to this point we have not specified the kernels 
/i 2 of the CS transformations 12101) or equivalently fixed 
the gauge for the statistical potentials a M and ai . The 
canonical choice 



/i,2(r — r') = ai,2 arctan 



y-y 



x ■ 



X' 



(13) 



with a\,2 = P, ±1, induces, through Fi, 2 (r — r') = 
lim £ ^ + a i,2% x (r — r')/(|r — r'| 2 + e), a symmetric 
gauge which is in conflict with the translational symme- 
try of the problem in the y direction. We are free, how- 
ever, to add to fi.2 functions 31,2 satisfying 31,2(1" — r') = 
31,2(1"' — r), corresponding to a regular gauge transforma- 
tion in the Lagrangian formalism. We choose 



3i, 2 (r-r) = lmai,a-= 

e ^ 0+ v(y-y) + £ 



x arctan 



V(y~y') 2 



(14) 



As a result (/ + .9)1,2(1" — r') = ^1,2(1 — e(y — y')] with 
e(x) — 0(x) — 1/2 where Q(x) is the step function. 
(/ + ,9)1,2 has a branch cut which we take to be along 
the positive x — x' axis and therefore V r (/ + 3)1,2 = 
[0, 2-kol\^(x — x')6(y — y')}. Accordingly, the statistical 
potentials [after the shift performed in obtaining i|10fl ] 
are in the Landau gauge 



a/(r) 



,(15) 



, - / dx' £ K u e{x - x')p,j{x', y) 
J J=i 

/m 
dx' £ K u e(x - x')pj(x', y)v y j(x, y) . 

(16) 



,7=1 



We will find it more convenient to obtain the static 
solutions and the edge modes in a slightly different gauge 
where we replace e(x — x') in (|15I16|I by 9(x — x'). Using 
this gauge we look for static solutions of the form 



Qi = Kiy — (j,i({m})t , pi = pi(x) 



(17) 



Taking for definiteness the confining potential to be an 
infinite wall situated at x < we set the density to zero 
on the edge. For the external vector potential we use 
the gauge A = [0, Bx] which gives velocity fields of the 
form Vj — and v\ — Vj{x). Assuming that <zj are 
also functions of x only and inserting the above ansatz 
into l|12afl - ljl2d(l we obtain a set of coupled equations for 
the various fields appearing in the problem. We solve 
the equations numerically for the cases of short-range 
interaction U = 27rA s (5(r — r') and Coulomb interaction 
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U = (27rA c )/|r — r'| with A c = e 2 / (2irelfrus c ) plus a con- 
stat positive neutralizing background. Representative so- 
lutions arc presented in Figs. ElandO! 



0.5 




FIG. 2: Mean-field solution for a v = 2/5 edge in equilib- 
rium K\ — k,2 — with short range interaction A s = f and 
Xij = 2nm/v5ij. The upper dotted and solid lines depict the 
densities pi — p2 and p = pi + p2 respectively, in units of po. 
The lower lines correspond to the current densities j\ — j% 
and j v = j\ + jf in units of poioj- 
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FIG. 3: Mean-field solution for a v — 2/3 edge out of equilib- 
rium K\ — 3, K2 = — 1 with Coulomb interaction— A c = 0.2 
and Xij = 2-Km/v5i.j. The upper dotted, dashed and solid 
lines depict the densities pi, p2 and p — pi + p2 respectively, 
in units of po- The lower lines correspond to the current 
densities j\, j'f and j y = j\ + jf in units of Pqlj c I. For in- 
creased asymmetry between the two condensates the strip of 
enhanced density near the edge becomes wider with filling fac- 
tor approaching 1, cf. Ref. 1331 The solution corresponding 
to the same parameters but when Ajj = 2-kKi.j is shown in 
the inset. 

The solutions to the mean-field equations form an 
m-dimensional manifold which is parameterized by the 



quantum numbers These quantum numbers are the 

conserved momenta along the edge of the different con- 
densates and may be thought of as their collective guid- 
ing center coordinates. Varying their values translates 
the condensates in the x direction. They also determine 
the amount of excess charge (per unit length in the y 
direction) carried by each of the condensates relative to 
a step-like constant density profile. To demonstrate this 
relation we integrate Eq. i|12a|) . noting that the ansatz 
(|17fl implies v V j (0) = kj and that the velocity vanishes in 
the bulk of the system. Consequently we obtain 

^2(K- 1 ) ij kj= / ( P/ (18) 
,i=\ J 

For a given set {«/} the values of {pLj} are determined 
by the requirement that far from the edge the density 
attains its bulk value v. To appreciate their physical 
meaning we consider the energy difference between two 
neighboring mean-field solutions differing by the number 
of particles in one of the condensates. For our ansatz 
such a variation is achieved by a suitable change of {ki} 
using H18fl . Utilizing the mean-field equations one finds 
in this case 

/m m 
^^^1=^^;, (19) 
1=1 1=1 

where J cPrSpi = SNj. Accordingly, /ij is interpreted 
as the chemical potential of the 7-th condensate. In the 
noninteracting limit we find numerically that the chemi- 
cal potential of equilibrium configurations (see below) de- 
creases and approaches hu> c /2 as the edge is placed away 
from the wall. This reflects the fact that under such con- 
ditions neighboring ground states may be related by the 
addition of charge to the rim of the electronic droplet. 
This result thus provides further evidence for the util- 
ity of the added contact interactions in the mean-field 
Hamiltonian. 

In case of short-range interaction the fields approach 
their bulk values exponentially in the distance from the 
wall2£ (cf. Fig. |2J). As a result we find from Eq. (|12d|) 
that for large x a/ n = — fij + v\ s + 1. Noting that in 
our gauge a/ (0) = we integrate Eq. Ijl2b|) to obtain 
the following relation connecting the chemical potentials 
{/x/} to the currents carried by the different condensates 



'£(K- 1 ) IJ » J --(l + v\ s )=I I . (20) 

,7=1 

For Coulomb interaction the asymptotic behavior of 
the fields at large distances from the edge is determined 
by the total charge carried by the system relative to 
the constant neutralizing background. From l|18|) we 
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find the linear charge density along the edge to be a = 
(v/m) 53r=i K i- Assuming that the length of the sample 
L is much larger than the length over which the density 
changes appreciably (~ cr/v) we can integrate the inter- 
action term in lfl2cil) to find (1/2tt) / U(r - r')[p(r') - 
v]d 2 r' ~ 2X c a\nL- 2X c f[p(x') - v\ ln(|x - x'\)dx' . The 
constant term is the leading contribution to the chem- 
ical potentials. It corresponds to the dominant part of 
the charging energy of the edge X c a 2 L\nL. The remain- 
ing contributions to the chemical potentials will be de- 
noted by flj — fij — 2X c ahiL. When u ^ we find at 



large distances pi ~ v/m+ 2X c a/x 



2X c a/x and 



ai ~ — fii + 1 — 2X c a\nx (with or without relative os- 
cillations between the condensates, cf. Fig 02). The slow 
decay of the current density results in a logarithmically 
diverging contribution to the integrated current from the 
wall to a point x. Defining Ii = Ii — 2X c a{v / m) lnx one 
obtains 



,7=1 



V 

m 



(21) 



If er = the fields relax to their bulk values more rapidly 



pi ~ vjm + 



and ai 



pi + 1 



Consequently the current is finite and a similar relation 
to lj2*T)) holds for {pi} and {//}. 

We remark that l|21l) leads to the appropriate quan- 
tized Hall conductance for a wide bar— {L ^> W 3> kI). 
Consider two edges in equilibrium with chemical poten- 
tials p\ t 2 and charge densities ±er. Summing over the 
edge channels we find 

I = h-h = v\pi ~P2 + 4A c( t MW/L)] = vVh*\\ ■ (22) 

The expression for the Hall voltage reflects the modifi- 
cation of the edge chemical potential due to the electro- 
static potential induced by the opposite edge. 

The particular solution to be realized out of the man- 
ifold of mean-field configurations is determined by the 
external constraints imposed on the system. Strictly 
speaking, since the ki are conserved quantities, the entire 
manifold corresponds to equilibrium states of the system. 
Taking, however, a more experimentally motivated point 
of view we will identify the equilibrium configurations as 
those which minimize the energy for a given total charge 
at the edge, as determined by the electrostatic potentials 
applied to the edge. Any perturbation which breaks the 
translation invariance in the y direction will cause the 
system to seek and eventually settle down in this solu- 
tion. A similar definition may be applied to the case 
where the edge is constrained to sustain a given total 
current instead of charge. 

As can be seen from Eq. I|19|) the symmetrical configu- 
rations (ki = k, pi = p) are extrema of the energy func- 
tional for a fixed total charge at the edge (J2i SNi = 0). 
From our numerical results we find that as long as this 
charge is not too large they indeed correspond to equi- 
librium solutions, see Figs. 01 and However, when the 
charge is increased the symmetrical solutions turn into 




FIG. 4: A map of the mean-field solutions for a v — 2/5 sys- 
tem with Coulomb interaction A c = 0.2 and Xu = 2nm/vSij. 
The solid curves correspond to equi-current families of solu- 
tions. The integrated current over a distance of 80/ from the 
edge is given in units of eu c . The dashed and dotted lines are 
the equilibrium (/ii = 112) and neutrality (ki — — K2) lines 
respectively. The reversal in the direction of the current with 
decreasing n signals a transition from skipping orbits along 
the wall to circular orbits on the rim of the Hall droplet. 



^ 






kJ 



FIG. 5: A map of the mean-field solutions for a v — 2/3 sys- 
tem with Coulomb interaction A c = 0.2 and Xu = 2nm/vSij. 
The notations are similar to Fig. 0] In the shaded region the 
mean-field solutions are unstable, as described in the text. 
The inset shows a similar map for Xu — 2nKu- Note the 
bifurcation of the equilibrium line ki = «2 upon entering the 
unstable region and the new asymmetric equilibrium states. 
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local energy maxima and new asymmetrical equilibrium 
states appear. This is particularly clear in the case of the 
v = 2/3 edge, see Fig. but the effect also exists for 
large k in the v — 2/5 edge. 

In the following section we will obtain the gapless edge 
excitations. As we have mentioned in the introduction 
they consist of a single magnetoplasmon mode and m — 1 
acoustic branches. For a v = 2/3 edge we find that the 
latter become soft inside a parameter region, denoted 
by the shaded area in Fig. whose extent depends on 
the particular model interactions (see also the inset in 
Fig. |BJ). Since in the acoustic mode the two condensates 
oscillate out of phase this softening signals a tendency 
of the edge to reconstruct by transferring charge from 
one condensate to the other, i.e. by deviating from the 
symmetrical solution along the direction of the neutrality 
line. Fig. confirms this conclusion. We note that in 
reality the edge reconstruction may involve population of 
additional Landau levels in the vicinity of the edge. Such 
processes, however, fall outside of the model studied here. 

When no constraints are imposed on the system it 
will choose the globally lowest available state. In the 
case of Coulomb interaction such a solution corresponds 
to a neutral configuration (k/ = 0) that minimizes the 
dominant charging energy (in realistic edges, for example 
in overgrown cleaved edge samples, additional potentials 
due to near-by charges may alter this conclusion 8 ). For 
short-range interaction the chemical potential of a system 
in equilibrium is a monotonically increasing function of 
k. Consequently energy is gained by placing the edge 
as far as possible from the confining wall. In order to 
fix the edge one has to introduce an additional confining 
potential. 

The equilibrium solutions compose only a small frac- 
tion of the mean-field manifold. The rest of the solu- 
tions describe possible realizations of the edge out of 
equilibrium. Deliberate action must be taken in or- 
der to place the system in one of these configurations. 
Experimentally, this may be achieved through selective 
population of edge channels employing gates over the 
edge. The relatively large equilibration lengths found 
in experiments^*^ (few /im) raise the possibility of ob- 
serving such states. 



IV. SYMMETRICAL EDGES 

A. Edge Excitations 

Having obtained the mean-field solutions we are now 
in a position to study small fluctuations about them. To 
this end we consider the RPA equations which are derived 
from Eqs. (|I2a|l - l|I2d|) by linearization around one of the 
mean-field configurations 
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FIG. 6: Excitations of a v — 2/5 edge in equilibrium 
ki = K2 = with Coulomb interaction A c = 0.2 and 
Xu = 2nm/vSij. The solid line depicts the density pro- 
file of the magnetoplasmon mode where 5p\ = 8p\- The 
dashed line corresponds to the density profile of the neu- 
tral mode Spi = —&p\ where the condensates oscillate out 
of phase. The upper inset presents the translation modes 
of this edge. The solid and dashed curves correspond to 
dpi/dni = dp2/dK2 and dpi/dK2 = dpijdnx respectively. 
The lower inset shows the velocity of the neutral mode for 
symmetric edges n\ = «2 = k with the same interaction. The 
solid and dashed lines correspond to a v = 2/5 and v = 2/3 
edges respectively. 



J2{K- 1 ) IJ e ij d i Sv j J = -Sp I , (23a) 
j=i 

m 

{K- 1 ) IJ (dtdiSdj - d t 5v l j - dMjo) = 

eiySpiv^ + e.jpjSvj , (23b) 

d t Spi = -V((5/7/t;/ + piSvi) , (23c) 
dtSOj - + P( PI , S PI ) + ^£ + Sai 

-^E A ^'-^ / U(v-r')5p{v')d 2 r' . 

(23d) 

Here P(pi,Spi) is the part of the linearized "quan- 
tum pressure" term, i.e., the second term on the right- 
hand side of (|12d|) containing x derivatives, and Sp — 

We can easily identify m independent solutions to the 
above equations. These are the small variations that con- 
nect the static solution, around which we linearized, and 
its neighboring configurations on the mean-field mani- 
fold. They are obtained from the static solution by taking 
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its derivatives with respect to the m conserved momenta 
Kj. Representative examples of these translation modes 
are shown in the inset of Fig. Evidently they cor- 
respond to a shift of the condensates in the x direction 
and hence to the addition and removal of charge from the 
vicinity of the edge. At present the amount of this charge 
is arbitrary since, due to the linearity of I|23a[l - lj23d|) . it 
is determined only up to a multiplicative constant. How- 
ever, we will show that the periodic boundary conditions 
impose on it discrete values in units of v /m. We there- 
fore interpret these modes as edge states of the m types of 
quasiparticles^ that can be generated by threading one 
of the condensates with an extra flux line [of the type 
associated with the CS transformation . 

The translation modes are the building blocks for the 
construction of the m gapless excitation branches sup- 
ported by the edge of the system. Considering the long- 
wavelength limit we find 



£■ a, k 

dpj' 



,7=1 



8k 



j 



,da Lo (x) 



IL 



dv\ (x) 



dnj 



u a 5 



L.J 



5Vj 



5 a 



y,a,k 



J=l 



dv V j(x) 



-i(ky—uj a t) 



(24) 



j=i Kj 



S9 C T 



k 



AJ(k)e 



-i{ky—Lj a t) 



where a = 1, ••• ,m is the branch index. The real and 
imaginary parts of l|24|) are, to first order in fc, inde- 
pendent solutions of the RPA equations. The expansion 
coefficients AJ(k) and the dispersion relations u a (k) are 
determined from the secular equation obtained by plug- 
ging (H3J) into (|23dJ) 



B 



E 

j=i 



[Bij 



U.J 



dn.j 



v 

2 A c — In 

m 



A a j(k) = 0, 



(25) 



7 is the Euler constant appearing in the long wavelength 
limit of the Coulomb interaction term (This and the fol- 
lowing results hold also in the case of short-range inter- 
action after replacing fl by \x and setting A c = 0). 

Due to the assumed symmetry of the mean-field solu- 
tion djli/dni = rji and dpi/ On j = 772 for all J ^ I. As 
a result the matrix B is symmetric and takes the form 



B = 



v 

2A r — In 



2e" 7 



C + (771 - 772) 1 . (26) 



Here again C is the pseudo-identity matrix. 

The eigenvalue problem l|25() is readily solved to give 
one edge-magnetoplasmon mode and m — 1 degenerate 
branches with acoustic dispersion^. The magnetoplas- 
mon corresponds to a coherent in-phase oscillation of the 
condensates with A\{k) = 2ir/Lv (cf. Fig. HJ). We have 
normalized the eigenvector such that the charge carried 
by its density profile is J d 2 rJ2T.j=iAj(k)(dpi/dKj) — 
27r (1 in dimensionful units) as can be easily verified using 
(|18p. Its dispersion relation is 



Vl 



k 



m 



1)172 + 2^A c ln 



2e^ 



(27) 



We will show below that this mode propagates in the 
positive y direction for all fractions. 

The rest of the modes are degenerate and have acoustic 
dispersion with velocity 



_ uj 2 

V2 = -r = m - m 



(28) 



They correspond to excitations in which the different 
condensates oscillate out of phase under the condition 
^2j—iAf(k) = and are therefore neutral (this perfect 
neutrality is not expected to hold in the case of a single- 
layer system). In the following we take them to be or- 
thogonal and normalized according to |A Q | = 2ir/L. 

The direction in which the neutral modes propagate 
depends on the filling factor. The reason for this stems 
from the relation between the velocity of the mode vi 
and the inverse compressibility of the edges of the con- 
densates. The latter can be viewed, according to Eq. 
(|18fl . as one-dimensional systems of length L and Nj — 
LY^j = \ {K 1 ) j j K J particles (where negative values of 
Nj correspond to minus the number of holes on the 
edge). Consider the inverse compressibility of the edge 
of the 7-th condensate xj 1 = {Nf /L)(Sfii /SNi) under 
the constraint of fixed total charge, i.e. for variations 
satisfying J2j^ K J = 0- Using Eq. I|18fl an d the form 
of K^ 1 one finds for the symmetric solutions xj 1 = 
±(Ni/L) 2 (r]i — 772). The stability requirement xi > 
then implies that sgn(v2) = ±1 for v = m/(pm± 1). 

We remark that for short range interaction the preced- 
ing analysis can be carried out to second order in k. To 
this order the RPA solution is the same as l|24|) but with 



" o ' fe multiplied by (1 + k). As a result 



opj , ovj and da 
one finds for the dispersion relations 

«i(fc) = [in + {m - l) m ]k{\ + k) , (29) 
wa(fc) = (m ~ V2)k(l + k) . (30) 



The inverse compressibility vanishes when the edge of 
the static solution is far away from the wall (i.e at large 
negative re) for which case u>2 ~ k 3 . A similar result for 
the case v = 1 was obtained in Ref. y3 
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B. The Tomonaga-Luttinger model 

In addition to the gapless edge modes a general ex- 
cited state of the system may include other types of 
excitations such as bulk magnetoplasmons, magnetoro- 
tons, quasiparticles and quasiholes. However, since all of 
them involve an energy gap the edge modes are expected 
to dominate the low energy physics. We will thus re- 
strict our consideration to cases in which the deviations 
of the various fields from their static ground state config- 
urations can be expanded solely in terms of the gapless 
excitations. Our plan is to plug these expansions into 
the Lagrangian l|l(J|) in order to obtain its low energy 
effective form and to use them to invert the CS trans- 
formations ll'Kil) . The result will be an expression for 
the original electronic operators in terms of the bosonic 
edge modes. For this purpose we are interested to use 
the kernels I|13I14I) and accordingly have to transform 
back to the gauge 1)151 llfijl which they induce. As a re- 
sult the gauge dependent quantities ai and Oi change, 
the latter to [«j — 1/2 J dxY^j = i ^uPj( x )]y ~ \pi ~~ 
1/2 J dxjyj = iKupjv V j(x)]t. The expansion of the de- 
viation fields is then 



5pi 

5v x ! 



5v^ 



Sa Io 



a— 1 k 

m 

pa( k ) Sv T a ' k i 

a— 1 k 

m 

EE P a (k)Svf a ' k , 



(31) 



a=l k 



E W°) 



]p a (k) 



K' + 5X>"(o)^ 



j=i 



Ok, 



Jo 



-iky 



1 rn 

a=l 



p a (o)A«(o) y + Y,p4We? 

k^O 



In (|31|) the time dependence has been shifted from the 
eigenmodes to the expansion coefficients p a (k) which are 
going to become upon quantization the creation and an- 
nihilation operators of the edge excitations. Similarly 
the coefficients p Q (0) will play the role of the number 
operators of the zero modes [the k — > limit of the 
real part of (|24(l ]. They, in turn, are related to the 
number operators pj(0) which measure the total excess 
charge at the edge relative to the background mean-field 
configuration due to the presence of translation modes 
(2irm/vL)(dpi /dnj). Explicitly 



Pi(0) = 



Lv v - 

27T771 2— 1 
a—1 



A?(0)p o (0) 



(32) 



In the following we take these operators to be time inde- 
pendent corresponding to fixed total charge at the edge. 



The effective low energy Lagrangian is obtained by 
substituting (|31|l into (|1(J|> and keeping terms to lowest 
order in k , 



' h ^TdtPi(k) Pl {-k) 



Lv ^— ' k 

fc>o 



2-nh m i 
T -j-J2^2-j-dtPa(k)p a (-k) 



L ^^k 

a=1 fe>0 

7T 



Lv 



[ m + (m~l)r, 2 }pi(0) Pl (0) 



tit, 

t j{m-m)^p a ^)p a {0) 



L 

yjj E ^i(k)pi(k)pi(-k) 

fc>o 



(33) 



L 



a=2 k>0 



where we have restored the units of dimensions. The 
first two terms of the Lagrangian constitute its symplec- 
tic part containing the information about the algebra sat- 
isfied by the density operators. Taking p a (k > 0) to be 
the coordinates we find that p a (k < 0) (up to constants) 
play the role of their conjugate momenta. Accordingly 
their commutation relations are given by {k > 0) 



[Pi(k),pi(-p)} = - 



kLv 



-8k P 



[p a (k),p (-p)} 
[pi(k),p a (-p)} 



2ir 
kL 

T-^—o a f3dkp 
Zn 



a,/3 = 2 > .-- 



= 



(34) 



The number operators p Q (0), or equivalently p/(0), are 
absent from the above algebra due to our assumption of 
constant edge charge. In case this assumption does not 
hold one needs to introduce additional operators, conju- 
gated to pi{0), that change this charge^. 

The effective low energy Hamiltonian is given by (mi- 
nus) the remaining part of the Lagrangian. Its first two 
terms correspond to the energy due to excess charge on 
the edge while the last two terms describe the energy of 
the gapless excitations. Using the commutation relations 
(|34|l one finds that p a {k) indeed create and destroy the 
eigenmodes (|24ll with frequencies uj a (k). Note that for 
the Hamiltonian to be bounded from below the velocity 
of the magnetoplasmon v\ should be positive while the 
neutral modes must propagate along the ±y direction 
(t>2 0) for the fractions v — m/(pm ±1), as argued 
previously and affirmed by our numerical calculations. 

To lowest order in k we find that in the gauge defined 
by (|15I16|) the bosonic operators are given by 4>i(r) = 

EBEE1 the CS 



y'piix) exp [i(9i + S6i)(y)]. Using Eqs. 
transformations J3 EJ) can be inverted to obtain from 
4>i(r) an approximate form for the electronic operator 
ipi(r). It contains the low energy edge components of 
the exact ipi(r). As a final step we project the two- 
dimensional operator on the edge by integrating over the 
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x direction with the weight d K pi(x)/ ' \J pi(x). The re- 
sulting expression for the one-dimensional operator [after 
omitting y-independent terms in A(r) and Aj(r)] is 

4>i(y) = L ~ 1/2 expi i( pi{y)} , 

m 

My) = Ky + J2 A i(o)p«(o)y (35) 

a=l 
m 

which constitutes the bosonization formula of the CTLL 
model. It is interesting that in the gauge we are using 
(pi is composed in equal parts of contributions from the 
CS phase A + Aj and from the phase S9j of the bosonic 
field. We also note, using H32I35|) . that the spectrum of 
pi(0) is constrained by the periodic conditions on ipj to 
take discrete values in multiples oiv/m. This observation 
conforms with our earlier identification of the translation 
modes as edge states of quasiparticles [see discussion fol- 
lowing pSaH23d|l ]. 

The above expression for ipi obeys fermionic anticom- 
mutation relations. It is a result of 

[<Pr{y)> Vj{y')] = (P± 5u)m sgn(y - y') , (36) 

which is easily obtained using the orthogonality rela- 
tion Ea l = 2 (^/) 2 = (27r/i) 2 (l - 1/m) (in order to have 
the proper anticommutation between different ipi one 
also needs to introduce Klein factors^.) Additionally 
Eqs. (|31I34|) can be utilized to derive [Sp(y), tpj(y')] = 

S(y-y')ipi(y')i wh ere Sp(y,t) = f dx YZLi S Pi( r > *)> in- 
dicating that creates a localized unit charge at the 
edge as it should. 

The bosonization formula (|35|) and the fact that the p a 
create eigenmodes of the system enable a straightforward 
evaluation of the edge propagator, which in the case of 
short range interaction is 

(4(y, t)^j(0, 0)) oc (y - Vl t)-i(y - v 2 t)-( 1 -^) ■ (37) 

The theory thus reproduces the power-law behavior 
I ~ V for tunneling into FQHE edges as predicted 
previously by Wen-i. One finds 7 = p + 1 for the 
v = m/(pm + 1) fractions and 7 = p + 1 — 2/m if 
v = m/{pm— 1). In the presence of Coulomb interactions 
the magnetoplasmon contribution to the propagator l|37|l 
introduces an additional correction ~ (lnt) - ^ at long 
times 2v\ c ij c t 1. 

So far we did not take into account the possibility of 
tunneling between the various condensates. Such pro- 
cesses may exist in multi-layer systems with small inter- 
layer separation (of the kind needed to establish correla- 
tions between the layers as assumed in our theory). They 
necessarily arise in single-layer FQHE systems from the 
terms connecting different composite fermions Landau 
levels as discussed following Eq. JSJ. In an m = 2 sym- 
metric edge Eq. (|35|l implies that tpi—i t 2 oc exp[i(y> c ± 



<p n )] where ip c and <p n are the magnetoplasmon and the 
neutral mode contributions to ipj respectively. Conse- 
quently a tunneling term ip\ip2 + h- c - modifies the neu- 
tral sector of the free edge theory (|33|l into a chiral sine- 
Gordon theory. The edge propagator for av = 2/(2p+ 1) 
system under such circumstances has been calculated by 
Naud et.aiSi and an algebraic decay, with the same ex- 
ponent as in Eq. (|37|l . was found to persist. 

V. ASYMMETRICAL EDGES 

When the edge configuration is asymmetric the matrix 
B appearing in the secular equation (|25|l is no longer sym- 
metric and the eigenvectors A" do not form in general 
an orthogonal set. More notably the latter do not corre- 
spond to eigenmodes of the effective Hamiltonian and an 
additional diagonalization procedure is needed. 

To demonstrate this point we proceed to evaluate the 
effective Lagrangian, as described in the previous section, 
continuing to use the normalization X)j=i = 2Trm/Lv 
in which the solutions carry a unit charge. The sym- 
plectic part yields the following algebra for the density 
operators 

kL 

[p a (k),p (-p)} = -— (M~ 1 ) a/3 5k P (k > 0) , 
/ T \ 2 m 

M a p=TP(pm±l)± I— ) J2 A i A i- ( 38 ) 

^ 71 ' 7=1 

Note that now A Q and therefore also M are (even) func- 
tions of k. For notational convenience we avoid an addi- 
tional k index on these quantities. The effective Hamil- 
tonian reads 

H = S £ d £-(K^) IL AtA^ Pam p m 

I..J.L 

a, (3=1 
a,/3=lfc>0 

It is now clear that the frequencies w a obtained from the 
secular equation are not the eigenfrequencies of H and p a 
do not create its eigenmodes. Nevertheless we can still 
express the tpi in terms of them. The result is similar 
to the one found in the symmetric case ll-'i3[l except that 
now k is replaced by kj. Using the identity proven in 
the appendix one can show that l|36|l and consequently 
the anti-commutation relations and the charge associated 
with ^1 remain unchanged with respect to the symmetric 
case. 

In the following we will derive the eigenmodes of l|39|l . 
consider their coupling to an external electric field and 
calculate the tunneling exponent into the edge. We will 
consider a generic fraction of the form v = m/(pm + 1) 
and restrict most of the discussion of the other fractions 
to v = 2/3. 
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A. v = m/(pm + 1) 

In order to make progress we note that the matrices 
M and K^ 1 are congruent 

M = (pm ± lfAK' 1 ! , (40) 

where A is defined in (|A.2|I . Accordingly they have the 
same signature. Since the spectrum of K~ x is composed 
of the eigenvalue v/m and m — 1 degenerate states with 
eigenvalues ±1 we conclude that M has m positive eigen- 
values for the v = m/(pm + l) fractions and a single pos- 
itive and to — 1 negative eigenvalues if v — m/ (pm — 1). 

Thus, there exists a (real nonorthogonal) matrix U 
such that UM^U = a where a is diagonal containing 
the signs of the eigenvalues of M~ l . Concentrating on 
the v = m/(pm + 1) case, where a = 1, we use it to 
define for k > the bosonic operators 




in terms of which 



Pa(-k) = 




V is an orthogonal matrix to be determined below. 
The hermitian conjugates of i|4H42|) are obtained using 
Pa(—k) — pj,(fc). Defining the real symmetric matrices 

N a p = ^M af3 (w a + up) , R = U- 1 NU- 1 , (43) 

the k part of H i|39|) takes the form 

a,/3=l Aj>0 

m 

Q=l fe>0 

provided we choose VRV = £1 with SI diagonal. Thus, 
once U and V have been calculated the eigenmodes of 
H are determined. In practice, however, we find them 
to be very close to the solutions of the secular equation 
()25J) . As a result, in the long wavelength limit and in 
the presence of Coulomb interactions the distinction be- 
tween a single magnetoplasmon mode and m — 1 modes 
with acoustic dispersion is preserved. The dispersion re- 
lations receive small, linear in k, corrections which lift 
the degeneracy between the acoustic modes. For short 
range interactions our numerical results indicate that in- 
creasing the amount of asymmetry tend to increase the 
velocity of the in-phase mode while decreasing the veloc- 
ity of the out-of-phase oscillations. 



Expressing tpj through a a k we can calculate the edge 
propagator. Assuming short range interaction we obtain 

m 

(4(y, t)M0, 0)) oc H (y - v a t)-^ , (45) 

a=l 

Using U~ 1 U~ 1 = M^ 1 and the result of the appendix 
one finds that the time dependence of the propagator is 
insensitive to the details of the edge configuration 

m 

71 = ^7/ =P + 1, (46) 

a=l 

in agreement with the universal tunneling exponent pre- 
dicted by Wenii. 
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FIG. 7: Mode couplings for neutral asymmetric configurations 
ki = — K2 = k of a v = 2/5 edge with short range interaction 
A a = 1 and Xu = 2-Km/v5i.j. The solid and dashed lines 
correspond to Q\ and Q\ respectively. 

For asymmetrical configurations the acoustic modes 
are no longer neutral. Consequently they couple to an 
external potential. The response of the edge to such a 
potential is given by the retarded density-density corre- 
lation function which is given by 

Z7T UJ — ilrvik) + 10 

a—1 v ' 

m 

The mode couplings obey the sum rule^ 

m m 

E^= E (^U = ^ (48) 

a—l a, 13=1 
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as shown in the appendix. Fig. depicts the way the 
couplings evolve as a v — 2/5 edge with short range in- 
teractions is taken out of equilibrium. Clearly the charge 
of the out-of-phase mode increases with the degree of 
asymmetry, at the expense of the in-phase mode. In 
the case of Coulomb interactions a similar trend exists 
but the couplings become fc-dependent. In particular the 
coupling of the acoustic mode decreases as 1/ ln(fc) in the 
long-wavelength limit (see also Fig. Illfl thus making it 
more difficult to detect. 
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FIG. 9: Phase velocities of edge modes in neutral non- 
equilibrium configurations ki = ~k,2 = k of a i/ = 2/3 edge 
with Coulomb interaction A c = 0.2 and A/j = 2nm/vSij. 
The solid lines depict, from top to bottom «i for k = 0, 1, 2, 3 
and 4. The dotted lines correspond in the same order to \v2\. 



with 



tanh20t = — - 



2R 



12 



Rll + i?22 



(51) 



FIG. 8: Edge excitations of a neutral v = 2/3 edge out of equi- 
librium ki — — k,2 = 2 with Coulomb interaction A c = 0.2 and 
\u = 2irm/vSij. They were calculated assuming a wavevec- 
tor kl — 10 -4 . The magnetoplasmon is shown in (a) where 
the dashed, dotted and solid lines present the density varia- 
tions Sp\,Sp2 and their sum, respectively. The corresponding 
quantities for the acoustic mode are shown in (b). 



B. 



2/3 



Here the edge supports two modes with a\ = 1 and 
02 = — 1. The introduction of the bosonic operators 



dak 



2tt\ 1/2 2 

X! UapPpi-Vak) , 
P=l 



kL J 



(49) 



Pa(-k) 



kL 
2^ 



1/2 r 



21: 



does not lead to a diagonalized H and an additional Bo- 
goliubov transformation 



h,2k = cosh.0fcai,2fc - sinhfl/caj 1 



k ■ 



(50) 



is needed to achieve the desired form. For an example of 
the resulting eigenmodes refer to Fig. [S] 

The eigenfrequencies are 



± (Jin - R 



22 



) + 7(i^7 



R 



22) 



-4Rj 2 



(52) 

While S7i 2 are both positive Eqs. I|49I5U|) imply that the 
two eigenmodes have velocities of opposite signs. The 
mode with acoustic dispersion continues to propagate 
with negative velocity i>2 = ~fl2/k. 

The velocities in the presence of Coulomb interaction 
are presented in Fig. The parameter k controls the 
width of the outer strip in which the density is larger 
than its bulk value v = 2/3 (cf. Fig. [3J). Varying its value 
introduces small 0(1/ In fc) corrections to the logarithmic 
k dependence of v\ and to the constancy of Vi- It does, 
however, modify the k independent parts of the velocities. 
Increasing the width decreases their values in accordance 
with the findings of Ref. |4(J. The numerical values found 
for v\,i are also quite close to those obtained in Ref. 0- 

The edge propagator in the presence of short range 
interaction exhibits a power-law behavior similar to l|45() 
but with non-universal exponents 
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a/3 



cosh26> fe + 2 (CT 



'02 



sinh 20i 



In Fig. HOI we present the tunneling exponents for neutral 
configurations with Ki — — K2 = k. Increasing n causes 
the leading edges of the two condensates to move away 
from the wall. When the edge nearest to the wall (that of 
the 1 — 2 condensate in this case) looses contact with the 
wall we find that 71.2 recover their equilibrium values. 
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FIG. 10: The tunneling exponents 71 (solid line) and 72 
(dotted line) for neutral non-equilibrium configurations n± — 
— K2 = c of a f = 2/3 edge with short range interaction 
X s — 0.2 and Xjj = 2-Km/i/Sjj. 

When tunneling occurs into the entire edge and not 
into a specific layer the long time (low energy) behavior 
is determined by the smaller exponent (72 in our case). 
Driving the system out of equilibrium has the effect of 
decreasing this exponent. While this tendency agrees 
with the experimental findings^^i of a smaller expo- 
nent than predicted by the CTLL theory in tunneling 
into FQHE edges, it does not provide a satisfactory ex- 
planations for the observed approximate \/v dependence 
of the tunneling exponent over a wide range of filling fac- 
tors including those of the form m/(pm +1). It should 
be noted here that in contrast to the symmetric case it is 
still unclear whether inter-condensate tunneling in asym- 
metric edges may affect the exponent. 

Finally, the couplings of the modes to an electric per- 
turbation along the edge, as deduced from the density- 
density correlator l|47[) are 

2 

Ql,2 - ( U ^)aia COSh6 k + ( U ^) a 2.1 S[nh ^ ■ ( 54 ) 

a=l 

They obey a similar sum rule to (|48ll 



(53) 

I 
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FIG. 11: The coupling of the acoustic mode for neutral asym- 
metric configurations ki = — K2 = ft of a v = 2/3 edge with 
Coulomb interaction A c = 0.2 and Xu = 2nm/vSij. The 
inset presents Q\ for similar configurations but in the pres- 
ence of short range interaction with X s = 0.2. According to 
15511 the coupling of the magnetoplasmon mode Q\ is given 
by2/3 + Qi 

2 

E ""Ql = v ■ (55) 

a=l 

The presence of Coulomb interactions makes the cou- 
plings acquire a 1/ In k dependence and reduces the cou- 
pling of the acoustic mode relative to the case of short 
range interactions (see Fig. 1 1 1 p . For the latter type of 
interactions we observe that once n is large enough such 
that the condensates are no longer in contact with the 
wall Q\ vanishes. The near neutrality of the counter- 
propagating mode makes its detection difficult especially 
in experiments^ where charge fluctuations are induced 
and detected over a wide range near the edge. The in- 
crease of the coupling of this mode to an external probe 
once the edge is taken out of equilibrium might help in its 
detection in similar experimental setups. Alternatively 
one should excite the two condensates selectively by care- 
fully placing a close top gate over part of the edge region. 



VI. CONCLUSION 

We have provided an explicit demonstration of the way 
in which the CTLL model for edges of multi-layer FQHE 
systems emerges as the RPA low-energy limit of the CS 
theory for such systems. The theory offers a unified de- 
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scription of edges both in and out of equilibrium and 
yields an edge spectrum composed of a single magneto- 
plasmon branch and m—1 acoustic modes together with 
their dependence on the underlying edge configuration. 
An attempt to carry the same treatment over to single- 
layer FQHE states is possible only if one avoids reference 
to the detailed structure of the composite fermion Lan- 
dau levels and neglects transitions between them. It is 
not clear to what extent such approximations are valid. 

As far as the question of the tunneling exponent is con- 
cerned, we find agreement with the universal results for 
the v — m/(pm + 1) states, obtained by previous studies. 
For v — m/(pm — 1) the minimal tunneling exponent de- 
pends on the structure of the edge and tends to decrease 
when the edge is driven away from equilibrium. While 
this behavior may account for some of the discrepancy 
between theory and experiment for v > 1/2 it is incon- 
sistent with the observed approximate \jv dependence 
of the exponent over a wide range 1 > v > 1/4 of filling 
factors containing both compressible and incompressible 
states. Moreover, additional processes, not taken into 
account here such as impurity scattering, may drive the 
system back to equilibrium and diminish the effect^. We 
hope that a careful study of the approximations needed 
in order to apply the theory to single-layer systems, as 
mentioned above, may help to settle this issue. 
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APPENDIX 

In this appendix we prove the identity 

(a=) E (M- 1 ) a ,A^ J =p±6 IJ , (Al) 



used in calculating the commutation relations and the 
propagator of ipj. To this end we define the matrices 

in terms of which we can write l|38|) as 



M = T — ^—rC±AA, (A.3) 
pm ± 1 

where C is the pseudo-identity matrix and tilde denotes 
transposition. Using the fact A a j = 1 we have 

AC = C and thus also A~ X C = C and CA^ 1 = C. As a 
result 



T = ±A~ 1 MA~ 1 = 1 ^—C . (A.4) 

pm ± 1 

Its inverse is readily found to be 

T- 1 = ±AKr 1 A=l±pC , (A.5) 

which is just Eq. 1A.1|I . Summing over / and J in Eq. 
(|A.ip we obtain 



m 

E (^ 1 U = ". ( A - 6 ) 

a, 0=1 



thus implying the sum rules (|48I55|I . 
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